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Large deflections of non-linear elastic rectangular plates 


By Lars WALLIN 


Introduction 


We consider a thin rectangular plate of a non-linear elastic material, hinged 
at the edges and loaded by a uniformly distributed transverse load p per unit 
area. The following investigation gives an approximate solution to the large, or 
membrane, deflection problem for the plate, neglecting the bending stresses. 

The plate has sides 2a and 26 and thickness h. We use a coordinate system 
with origin at the midpoint and the xy-plane in the middle plane of the plate. 
The z-axis is parallel to the side of length 2 a, and the z-axis is perpendicular 
to the plate and positive in a direction opposite to the load p. The displacement 
has the components w, v and w in the coordinate directions. The plate has a 
stress-strain relation having the same analytical expression as Norton’s creep 
law, namely 


BiG 4 


or generalised to three dimensions [1] 


with 


with 


The stress invariant o, and the strain invariant e, are 
ZS 
Oe = 5 S55 Sif 
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A plane stress distribution is assumed and the stress tensor, the stress deviation 
and the stress invariant then are 
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The plate is assumed incompressible (v=1/2) so we have e,+é,+6,=0. With 
the notation ¢,,=1/2y1, where y,, is the usual shearing strain, and neglecting 
Yr, and y,,, the strain tensor and the strain invariant are 
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Only deflections larger than the plate thickness are considered, and the bending 
stresses can be neglected when compared with the membrane stresses. The rela- 
tionship between strain and displacement then is (reference [2], p, 304.) 
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Energy method solution 


For the rectangular membrane the equations of equilibrium together with a 
relationship between stress and displacement are very complicated, but it is 
possible to determine an approximate solution to the membrane problem using 
an energy method. A solution is assumed with a number of arbitrary constants, 
which are then determined so that the potential energy is a minimum. The 
general expression for the potential energy for the transverse loaded plate is 


U=({ff [ox de, + oy dey + 2ondey|dadydz+ |{ pwdady—f ~(P6+MO6)ds (6) 


The first term is the strain energy, the second term is the work of outer forces 
and the third term is the work of the reactions. The third term is zero for a 
simply supported or a clamped plate, and is omitted. The strain energy per unit 
volume can be transformed by (2), (3) and (4) and integrated. The z-integration 
in (6) can also be carried out, and then the potential energy becomes 
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U=K IJ e*dxdy+{{ pwdedy (7) 
where 
K- n Ah 
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The integration is carried out over the plate area. 
For the simply supported plate we assume the solution 
w= A (x? —a?) (y? — Bb?) 
u = Ba (x* — a”) (y® — b*) (8) 
v= Oy (x? —a?) (y®— 8) 


eA [ac (y* b?)? 4 | yf (a? a) 
1 


pm ie (et =a) Ue O(a) (ye 0) 
+38 A? O[y? (3 y? — 6?) (2? —a?)? + £2? (5 y? — B?) (y? — B?)?] (a® — a?) + 
+4 B?[ (3.27 — a?) (y? — B*)? + 2? y? (2? — a*)?] + 
Fee ley —0) (=a) 4e yyy 
+ebCl((3a —a) (sy —b)+2a7y"] (2° —a) (yy? —0") (9) 


It is not possible to integrate (7) for x+1. In formulating the minimum 
conditions 


(10) 
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we have to change the order of differentiation and integration, and after some 
rearrangement these conditions can be written 
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From a dimensional analysis of (8), (9) and (11) we find that it is possible to 
express the constants A, B and C as 


os n+2 
1 ( palk 
A a NG ; 
a h (12) 
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where the coefficients «, f and y depend only on the ratio b/a=/. We define 
a polynomial H as 


H=4F,+4p@F,+4y@F,+f'a'F,+7'a'F,+Bya'F, (13) 

so that 
ee=gA°H (14) 
The polynomials F are obtained from (9), (13) and (14). With (12), (13) and (14) 


we can rewrite (11) and get 
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The three unknowns «, 6 and y can now be calculated from (15) with a digital 
computer. The results of such calculations are given in table 1. The maximum 
deflection of the plate is 


“== ada 22) ‘ (16) 


The linear case 
In the special case n=1 (x=1) we can integrate (7) exactly and get 
U=K(P,A*+ P, A* B+ P, A?0+ P, B+ P,0?+ P,BC) + pa®b? A 
where P, to P, are 
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Table 1. 
bla=A n B y oO wae 
1 i 0.442 0.442 0.694 0.694 
3 0.446 0.446 0.494 0.494 
5 0.449 0.449 0.426 0.426 
7 0.451 0.451 0.392 0.392 
10 0.452 0.452 0.365 0.365 
co 0.456 0.456 0.293 0.293 
1.5 1 VL72 0.280 0.375 0.844 
3 1.026 0.419 0.280 0.630 
5) 0.963 0.463 0.246 0.554 
7 0.928 0.485 0.229 0.516 
10 0.896 0.504 0.215 0.485 
co 0.801 0.553 0.178 0.401 
2: 1 2.195 0.0863 0.226 0.904 
3 1.856 0.370 0.173 0.691 
5 1.684 0.450 0.153 0.613 
7 1.578 0.490 0.144 0.574 
10 1.481 0.523 0.136 0.542 
co 1ST 0.613 0.114 0.457 
3 ul bLO2 — 0.322 0.105 0.944 
co 1.858 0.684 0.0559 0.503 
co 1 0.96 
56 - 128 
= 2 ) h2 525 
ig. 51-25-63 126°) a’b 
56-128 
12. = b? — 12a7)a°b° 
3° 21-25-63 ( ) 
7°128 ee 
Eyam + 84.5”) a? b° 
7 ° 128 9 2 875 
= tty APA / 5 
ee 51-95-63?” + 84a") a’ b 
1A ie 128" 2: 
eee 
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The minimum conditions (10) then give 
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For the square elastic plate with »=4 we obtain from table 1 or from (17) 
ao =0.694. For y=1/4 we can find the value «=0.802 in [2] p. 346. From (17) 
we can determine the maximum deflection for a very long plate (b>qa). In this 
case the quantity «A in (16) will be equal to 0.96. 


Simplified solution 


For the non-linear case («+-1) we have to choose a simpler solution than (8) 
to be able to evaluate it without a computer. We obtain such a solution by 
putting B=C=O in (8), and the energy expression then becomes 


U=4K (fy A*- I+ part? A 
with 


ab 
T=[ [Wyre ty (ea) dedy 
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The integral J must be evaluated numerically, but this can easily be done by 
hand-calculation. For the constant A we get 


n 


ie n+1/(3\*paebVk igen a a 
a n \16) 9xhI ak ps 


In the case of a square plate this simpler solution gives a maximum deflection 
which is about 10% less than the one obtained by the complete solution (8). 


Discussion of the accuracy 


To get an impression of the accuracy of the method used above, we can 
apply it to the corresponding problem for the circular plate, where an exact 
solution is determined by Odqvist [3]. For the circular plate we assume 


where r is the radius vector, a the radius of the plate and uw and w the radial 
and transverse displacement. For n= 1, which case is easy to calculate, the energy 
method gives a maximum deflection, which is 1 or 2% too large. A simpler 
solution with B=0 gives a maximum deflection which is between 10 and 14% 
too small. It is reasonable to assume that the deviation from the exact solution 
is somewhat larger for a square or rectangular plate than for the circular plate. 
The accuracy is also less for longer plates. Other examples of the accuracy of 
the energy method can be obtained from small deflection problems with n=1, 
for instance from the problem of a uniformly loaded clamped rectangular plate, 
where a comparable energy method [4] gives about 5% error for b/a=1, 12°% 


error for b/a=2 and about 25% error, when b/a becomes very large. All these 
values are above the exact value. 
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